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Abstract 

Simulations have shown that as two metal surfaces approach 
each other, the surface layers can avalanche together when the rigid 
interfacial spacing falls below a critical distance. This is accompa- 
nied by a discontinuous decrease in the adhesive energy. Here we 
present an examination of this phenomenon for the bcc metals Fe 
and W using the Equivalent Crystal Theory. In order to identify the 
circumstances under which avalanche might be inhibited, the effect 
of loss of registry between the two surfaces is investigated in detail. 
The avalanche is inhibited when the two surfaces are sufficiently far 
out of registry and when only a few layers near the surface are al- 
lowed to relax. As the relaxing slabs get thicker a sharp avalanche 
reappears. However, as the loss of registry increases the energy 
released in the avalanche decreases. 

PACS Numbers: 68.35.-p, 46.30. Pa, 81.40.Pq 



I. INTRODUCTION 


Recent simulation studies using different semi-empirical methods 1-6 suggest 
that there are conditions under which solid surfaces will jump across and close an 
interfacial gap even when the initial interfacial separation is significantly larger 
than the bulk interplanar spacing. The first suggestion that solid surfaces could 
jump together was made by Pethica and Sutton. 1 Primarily interested in the 
interaction of a metallic tip with a flat surface in the scanning tunneling 7 and 
atomic force 8 microscopes (STM/AFM), they drew their conclusions from studies 
using either a Lennard-Jones pair potential or continuum elasticity theory. Both 
of these approaches have serious limitations of which the authors were certainly 
aware. 1 Subsequent investigations 2-6 of the stability of adhering flat surfaces have 
found quantitative evidence for an avalanche effect. 

Smith et al. 2 studied a (100) interface between two nickel half crystals using 
the Equivalent Crystal Theory (ECT) 9,10 which is based on the universal binding 
energy relation (UBER). 11 Initially the atoms in each half crystal were fixed rigidly 
at bulk interatomic positions and the adhesive energy was computed as a function 
of dR, the rigid in ter facial spacing. Next, at each value of dp_ a gradually increasing 
number of planes parallel to the interface were allowed to relax in order to minimize 
the total energy. Allowing only the surface plane on each slab to relax, it was found 
that for large values of da , consistent with relaxation at free metal surfaces, 10,12-14 
the surface layer on either half crystal relaxed inward reducing the separation 
between itself and the first subsurface layer. There was a significant energy barrier 
preventing the surface layers from jumping across the interfacial gap. As the value 
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of d,R was decreased, the height of this energy barrier decreased until at a certain 
critical value of dfi it disappeared. At that point it became energetically favorable 
for the two surface layers to jump across the interfacial gap and come together. In 
other words, when the two surfaces were pushed to within some critical distance, 
djP*, of each other, the surface atomic layers could not be held apart. This jumping 
across to close the interfacial gap was rather sudden and was accompanied by a 
sharp, discontinuous drop in adhesive energy. The large residual elastic strain, 
now distributed over three interatomic spacings rather than being concentrated 
at the interface, was gradually relieved as the distance betw'een the substrates was 
decreased farther. As a result, the post-avalanche variation of the adhesive energy 
with <1r was more nearly parabolic. This was quite different from the behavior of 
the energy in the case of rigid adhesion which followed the UBER 15 and is well 
described by the Rydberg function: 

E(d) = AEE*(d*), (4a) 

E*(d*) = — (1 + d*)exp(—d*), (46) 

d* = (d-d m )/l, (4c) 

Upon relaxing more than one atomic layers it was discovered 2 that as the number 
of relaxed atomic layers, n, increased, the critical value of dR for avalanche to 
occur also increased. Smith et al 2 also argued that as n — ► oo, the critical 
distance d^ 11 ~ Inn and the discontinuity in the energy approaches the surface 
energy which is of the order of 1 eV per surface atom. 

In their first report Smith et al. 2 indicated some circumstances under which 
avalanche at solid surfaces may be inhibited. Among other things, they suggested 
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that a lack of registry across the interface may prevent, or at least severely inhibit 
the avalanche. In a subsequent zero-temperature Monte Carlo (MC) investigation 
of the avalanche phenomenon in Ni Good et a/. 6 reported that in the case of Ni 
(001) a total loss of registry - so that atoms on either side of the interface come 
down on top of each other - avalanche is severely inhibited in that there is no sharp 
drop in the energy. However, there was a significant, if rounded, avalanche-like 
change in the interlayer spacings. That study 6 allowed no more than three surface 
layers on either side to relax and the authors noted that as the number of relaxing 
layers increased the changes in the interlayer spacings appeared to get sharper. 

In this paper we present a study of avalanche in bcc metals. We have investi- 
gated avalanche at the (001) and (110) interfaces of Fe and W. Although W (001) 
is known to reconstruct, we have ignored reconstruction effects for the present. 
We also present here a more detailed study of the effect of registry on avalanche 
at the (110) interface of Fe. In the next section we present a brief review of the 
simulation procedure while in Section III we present the results for avalanche at 
the (001) and (110) interfaces in registry for both, Fe and W and the results for 
the out-of- registry (110) interfaces of Fe. Finally, we summarise the study and 
discuss directions of future work in Section IV. 

II. SIMULATION PROCEDURE 

In the simulations that we present here we have used the Equivalent Crystal 
Theory 9,10 (ECT) to compute the energies of the system. This method, based on 
the Universal Binding Energy Relation, 11 has been described in detail elsewhere. 
Nevertheless we present here a very brief review of the essentials. The ECT, as 
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generally implemented, expresses the energy of a collection of atoms as a sum 
over individual atomic contributions. Each atomic contribution comprises four 
different terms. The first of these terms depends essentially on the local density 
in the immediate neighborhood of the atom in question and is generally the largest 
single contribution to the surface or interface energy. The second term accounts 
for local deviations in symmetry away from that of the ground-state crystal and 
local variations in nearest-neighbor distances. The other two terms depend on 
changes in bond angles and account for shear-like distortions. The ECT has 
been shown to give accurate surface energies and surface relaxations for a variety 
of materials. 8,10,12,13 In particular, it has been shown 13 that the last two, bond- 
angle-dependent terms contribute little to the relaxation energies of metal surfaces. 
Hence, in this study we have neglected the last two terms of the ECT energy 
expression. 

This study has been restricted to planar relaxations only, t. e., each atomic 
layer parallel to the interface is assumed to retain its planar structure as the system 
relaxes - no in-plane reconstruction or buckling of atomic planes is permitted. In 
order to relax the system the energy was minimized with respect to the various 
interlayer spacings using a Monte Carlo procedure. The energy of the system - 
two semi-infinite solids with parallel surfaces either in or out of registry - was 
calculated at each value of the rigid interfacial separation, dn, as a function of the 
interplanar spacings in either half solid. The variations in interlayer spacings were 
assumed to be symmetric about the interface. That is, the distance between layers 
i and t + 1, di,i+i, had the same value in the two half solids. Hence, when allowing 
n layers to relax one needed to be concerned with only n independent variables 
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and to compute the energy of n 4- 2 layers in only one of the semi-infinite slabs. 
As a result the energies presented here are half of the total system energies of 
adhesion and the depth of the energy well is the surface energy of the appropriate 
metal surface. 

It should be noted that these calculations do not make use of the complete 
Metropolis Monte Carlo algorithm 16 in that the temperature of the system is as- 
sumed to be zero and there is no provision made here for escape from a local 
energy minimum. The energy surfaces are rather smooth and the general loca- 
tions of the energy minima are known. As discussed earlier, for a given d.R there 
is either one minimum (for those values of cLr for which avalanche does occur) 
or two minima (for those values of cLr where avalanche does not occur). One 
of the minima occurs when all relaxing layers are closer to each other than to 
their repsective bulk slabs. In the non-avalanche case, there is a second (local) 
minimum which occurs when ail relaxing layers are close to their respective bulk 
slabs. Avalanche occurs when the energy barrier between this minimum and the 
global one with the relaxing layers closer to each other disappears and the local 
minimum gives way to an inflection point. It is precisely the transition from a 
two-minimum configuration to a single-minimum one (as dn is decreased) and the 
associated transition of the system from a local energy minimum to a global en- 
ergy minimum that we wish to investigate. This is why we have not implemented 
the full Metropolis MC algorithm. However, we do plan to implement the full 
Metropolis MC scheme soon - allowing individual atoms their freedom in a larger 
computational cell instead of requiring a planar movement - in order to investigate 
the effect of temperature on the phenomenon of avalanche. 
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III. RESULTS 


A. Interfaces in registry 

The variation of relaxed binding energy with initial rigid interfacial spacing, 
d,R, for the case with the approaching surfaces in perfect registry, is shown in 
Figs. 1 and 2 for the (001) and (110) interfaces of Fe and W. The sudden drop in 
the energy upon avalanche is clearly evident. 

The physical transformation associated with this precipitous drop in energy 
is seen clearly in Figs. 3-5. Figure 3 shows - for different interfaces (in perfect 
registry) and different numbers of relaxing layers as indicated - plots of the final 
separation between the relaxed surface layers on either side of the interface, also 
measured with respect to the equilibrium interplanar spacing, as a function of 
cLr. At large values of dR, when the surfaces are essentially isolated and not 
interacting, the separation between surface layers is actually slightly greater than 
dft. This implies that the surface layer relaxes in toward the bulk - a phenomenon 
known from experiment as well as theory. 10,12-14 As dR decreases, this separation 
decreases slightly, indicating a gradual reduction in the inward relaxation of the 
surface layer and indeed even a slight outward relaxation as the interaction across 
the interface increases in strength. At a certain “critical” value of dR, however, 
there is a precipitous drop in the separation between surface layers, as they leave 
their respective slabs and come together in the center. 

In Figs. 4 and 5 we have plotted the values of the different interplanar sep- 
arations, di,i+i against dR for some of the cases shown in Fig. 3. Here we see 
essentially a repeat of the phenomenon seen in Fig. 3. Although the effect is hard 
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to detect on the scale of these figures, for large values of dp, d ti , + i is greater 
or smaller than the equilibrium (bulk) value depending on whether t is even or 
odd. This is a reflection of the well-known damped oscillatory nature of surface 
relaxation. 10,12-14 Again, at a critical value of dp. there is a dramatic change, this 
time an increase in <2i,i+i indicating that the relaxed layers have separated from 
their respective slabs and form an elastically strained zone between the two bulk 
slabs. When more than one layer is relaxed, it is curious that while the changes in 
almost ail the dt,»+i are sharp increases followed by an almost linear decrease, the 
behavior of dip and dn,n+i show a slower rounded increase after the initial steep 
rise before they too settle into a near-linear decrease to the equilibrium value. 

As one of the objectives of this study was to investigate the behavior of d < ^' t , 
the critical value of dp for avalanche to occur, with the thickness of the relaxing 
slab, we have plotted in Fig. 6 the variation of df^ 1 with n, the number of relaxing 
layers. As the exact point of avalanche in our simulations is somewhat dependent 
on the history of the simulation, the values of d^ tt are impossible to determine 
precisely. Hence, we have attempted to estimate an uncertainty in the values of 
d ft* 1 which is reflected in the error bars on the points. Anticipating a logarithmic 
variation of df^ 1 with n, as discussed in the Introduction, we have also plotted a 
least squares fit to a logarithmic function to the data. The actual values of the 
parameters and the regression are indicated on the graphs. It might be argued that 
the quality of fit is not excellent. We believe this is because the basis for expecting 
a logarithmic variation - arguments involving continuum elasticity theory and the 
primarily exponential tail of the force vs. dp curve in the case of rigid adhesion 
- is sorely tested in avalanche involving only a handful of atomic layers near the 
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interface. As the thickness of the relaxing slab increases, one would expect to find 
a better fit to a logarithmic function. 

B. Interfaces out of registry 

So far we have presented results for avalanche at two different interfaces, (001) 
and (110), of two different metals, Fe and W, when the two approaching surfaces 
are in perfect registry with each other so that when the value of d,R is reduced to 
zero the interface disappears and the two halves join to make a complete, infinite, 
solid. We have seen that the results are very similar for the two metals and 
for the two interfaces. In this section we present results for cases when the two 
approaching surfaces are not in registry, i. e., starting from the case of perfect 
registry, one half solid has been shifted relative to the other by small amounts in 
a direction parallel to the interface. 

In a previous study 6 we had presented preliminary results for avalanche in 
the case of out-of-registry approach for the (001) interface of Ni where the shift 
had been by half the distance between neighboring surface atoms so that atoms 
on either surface layer came down on top of one another as da was reduced to 
zero. In that paper we had also allowed up to three layers to relax. We had 
reported that the ’’total” loss of registry strongly inhibited avalanche. However, 
we had pointed out indications that the transitions, which were rounded rather 
than sharp for total loss of registry, grew sharper as the number, n, of relaxing 
layers increased. This raised the possibility that if n increased to macroscopic 
numbers avalanche would no longer be inhibited. In this section we investigate 
this phenomenon for the Fe(110) interface. We do not present results for out-of- 
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registry W(110) separately because, as we have already seen for the in-registry 
case in the previous section, they are essentially similax to those for iron. 

Figures 7 and 8 show plots of relaxed interfacial spacing and relaxed adhesive 
energy, respectively, against dp for out-of-registry Fe (110) interfaces with one and 
five layers relaxing. Results are shown for different amounts of shift, as labeled, in 
the ( 1 10) direction. We can see from Fig. 8 that for the cases of shifts intermediate 
between perfect registry and total loss of registry there appears to be no sharp drop 
in energy associated with avalanche when only one surface layer is allowed to relax. 
However, upon allowing more layers to relax we do see a sharp discontinuity in 
the energy associated with avalanche. Figure 7 shows the associated change in the 
relaxed interfacial separation. There is a sharp discontinuity in the rate of change 
of the interfacial separation with dp. associated with the discontinuity in adhesive 
energy. This is made clearer in Figs. 9 and 10 where we have plotted the values 
of the interlayer spacings as functions of dp for the situations depicted in Figs. 7 
and 8. Figure 9 shows results for one layer relaxing while Fig. 10 shows results for 
the case of 5 layers relaxing. We see that even with only one layer relaxing there 
is a sharp change in the behavior of the interlayer spacing at a critical value of dj? 
although the subsequent variation of d \ t 2 is quite different for the out-of- registry 
cases from that seen in the case of perfect registry. This indicates that there is 
avalanche of sorts occuring even with one layer relaxing. From the small size of 
the drop in energy seen with five relaxing layers, we conclude that there might be 
an energy drop even in the one-layer case but that it is too small to be observed 
in our computations. 
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IV. DISCUSSION AND SUMMARY 


The phenomenon of avalanche can be viewed as a competition between the 
attractive interaction of a surface layer with the layer(s) directly below in the 
same slab, and with the surface layer(s) across the interface in the other slab. At 
large values of d/j, the surface layer only interacts with other layers in the same 
slab, and thus it moves in toward the bulk of the slab. As cLr decreases, however, 
the surface layer begins to experience an attractive interaction with the layers on 
the other side of the interface. Eventually there comes a point at which forming 
bonds across the interface lowers the energy sufficiently to compensate for the 
increase in energy caused by stretching the bonds between the surface layer and 
its neighboring bulk layer. At this point avalanche occurs. 

The suppression of the effect in the totally out-of-registry case for thinner 
slabs, at least, can be understood by looking at the surface geometries. In the 
in-registry case, each (110) surface atom is bonded to two atoms in the neigh- 
boring bulk layer, four atoms in the same layer, and (for small enough cLr) two 
atoms across the interface. Thus the surface atom will begin to interact with 
the two atoms across the interface at the same time, and the interactions will 
be equal in magnitude as long as perfect registry is maintained. On the other 
hand, upon total loss of registry, a surface atom is bonded to two atoms in the 
neighboring bulk layer, and four in the same layer, all within the same slab. In 
the other slab, however, there is only one “nearest-neighbor” atom directly across 
the interface. As the slabs approach each other, the interaction across the in- 
terface is substantially weaker (relative to interactions with atoms in the same 
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slab) than in the in-registry case. There are two other atoms in the surface plane 
across the interface which might be called “secondary” near neighbors rather than 
second-nearest neighbors. The interactions with these, however, will be strongly 
screened, especially at small values of cLr. Consequently, these cannot compensate 
for the decreased coordination across the interface. This is what inhibits the sharp 
avalanche in the out-of-registry case. 

As the number, n, of relaxing layers increases, however, the deeper layers play a 
significant role and a significant amount of energy can be gained by relaxing these 
outward. Now, when the surface layer moves out in order to “gain” coordination 
the increase in energy caused by stretching the bonds with the substrate can be 
minimized by distributing the “strain” over severed interlayer bonds. This is what 
causes a return of avalanche when more layers are allowed to relax. This leads 
us to believe that for experimentally relevant thicknesses loss of registry will not 
totally suppress avalanche. However, as is evident from a comparison of Figs. 1 
and 2 with Fig. 7 or Figs. 4 and 5 with Figs. 8 and 9, the value of does 
decrease as we move farther out of registry. Hence, loss of registry does appear to 
inhibit avalanche without suppressing it totally. 

In summary, we have carried out a Monte Carlo investigation of avalanche in 
Fe and W. In agreement with previous work, 1-6 we find that an adhesive avalanche 
effect may be observed when two metal surfaces are brought into proximity. When 
the interfacial gap between these surfaces falls below a critical value (which de- 
pends on the number of layers relaxed) , the relaxing surface layers will move away 
from their respective bulk substrates and form a strained zone between them. 
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We have extended these previous results by considering in detail the effects 
of registry on the avalanche phenomenon. We find that the avalanche transition 
is sharp when the two surfaces are in exact registry. When the surfaces are out 
of registry with each other the transition is no longer as pronounced but that 
avalanche occurs nonetheless if several layers on either side are allowed to relax. 
We further find that the critical distance for avalanche, d < p t , and the energy 
released upon avalanche both decrease as the loss of registry increases. 

It must be pointed out here that no lateral movement parallel to the interface 
was permitted in these simulations, It is possible that when several layers are 
allowed to relax both perpendicular and parallel to the interface they may move 
laterally in order to restore registry and thus facilitate the avalanche process. 
This would introduce an additional shear-type strain in the relaxing zone but 
this increase in energy may be more than compensated for by the energy released 
upon avalanche. This is an aspect of avalanche into which we are currently looking. 
We are also currently in the process of incorporating non-zero temperatures into 
atomistic simulations on larger unit cells as we have already mentioned in the 
Introduction. 
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Figure Captions 


Figure 1. Plots of relaxed energy versus rigid interfacial separation, dp, for Fe 
(a) (001), (b) (110) interfaces. 

Figure 2. Plots of relaxed energy versus rigid interfacial separation, dp, for W 
(a) (001), (b) (110) interfaces. 

Figure 3. Plots of relaxed interfacial separation versus dp for (a) Fe(001), (b) 
Fe(110), (c) W(001), and (d) W(110) interfaces. 

Figure 4. Plots of interplanar spacings, di,»+i for Fe: (a) (001) interface with 1 
layer relaxing, (b) (001) interface with 5 layers elaxing, (c) (001) inter- 
face with 7 layers relaxing, (d) (110) interface with 5 layers relaxing. 

Figure 5. Plots of interplanar spacings, di,i+i, for W as in Fig. 4. 

Figure 6. Values of d ^ nt vs. n, the number of relaxing layers, plotted for (a) 
Fe(001) and Fe(110), and (b) W(001) and W(110) as labeled. The 
solid lines show least-squares fits to a logarithmic function. The fitted 
functions and the vaues of the correlation are shown on each plot. 

Figure 7. Plots of relaxed interfacial separation versus dp for Fe(110) interfaces 
out of registry. One surface has been shifted relative to the other in the 
(hkl) direction by varying amounts - in units of dhkl, the equilibrium 
interplanar distance - as indicated on each of the graphs. Results 
shown are for cases with 1 and 5 surface layers relaxing. 

Figure 8. Relaxed energy vs. dp plotted for Fe(110) interfaces out of registry. 

One surface has been shifted relative to the other by varying amounts 
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- in units of dhki , the equilibrium interplanar distance - as indicated 
by the vector b on each of the graphs. In each graph results are shown 
for cases with l and 5 layers relaxing as indicated in the legend. 

Figure 9. Plots of interplanar spacings, d,,i+i for Fe(110) interfaces out of reg- 
istry. One surface has been shifted relative to the other in the (hkl) 
direction by varying amounts - in units of lattice constant - as indi- 
cated on each of the graphs. Results shown axe for cases with 1 surface 
layer relaxing. 

Figure 10. Plots of interplanar spacings, dt.i+i for Fe(110) interfaces out of reg- 
istry. One surface has been shifted relative to the other in the (hkl) 
direction by varying amounts - in units of dhki, the equilibrium inter- 
planar distance - as indicated on each of the graphs. Results shown 
are for cases with 5 surface layers relaxing. 
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